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Abstract 

We present constraints on the CKM parameter 7 = arg V*i^ formed within the framework of 
SU(3) symmetry and based on charmless hadronic B decays to K*^Tr^ and other pseudoscalar- 
vector final states. For strong phases of 0(10°), our analysis weakly favors cos 7 < 0. We also 
estimate that a determination of 7 with an experimental uncertainty of less than 10° can be attained 
with an order-of-magnitude improvement in the precision of the experimental inputs, but SU(3) 
symmetry breaking could introduce corrections approaching the size of the current experimental 
uncertainties. 
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In the Standard Model, the Cabibbo-Kobayashi-Maskawa (CKM) quark mixing matrix [l| 
gives rise to CP-violating phenomena through its single complex phase. This phase can be 
probed experimentally by measuring decay rates and CP asymmetries for charmless hadronic 
B decays that receive contributions from amplitudes with differing weak phases. In the 
flavor SU(3) decomposition of the amplitudes in pseudoscalar- vector {PV) final states 0, 
the h — i> uus transition B — >■ K*'^'k^ is dominated by two amplitudes, color-allowed tree and 
gluonic penguin, which interfere with a weak phase vr — 7, where 7 = arg V*fj, and with an 
unknown strong phase 5. We extract cos 7 in two ways, employing Monte Carlo simulation 
to propagate experimental uncertainties and ratios of meson decay constants to account for 
SU(3) symmetry breaking. The first uses only B —>■ K*^7[^, and the second adds information 
from B — i> 0'^* . In both cases, the magnitudes of the penguin and tree amplitudes must be 
known, and we estimate these from CKM unitarity and measured branching fractions for 
other B PV decays. 

An alternative method of constraining 7 makes use of observables pertaining to 6 — c 
transitions and to mixing in the neutral B and K systems, with a resultant 95% confidence 
level (C.L.) allowed interval for 7 of [38°, 80°]. In contrast, the analysis presented in this 
Letter uses rare, charmless b —>■ u,d, s transitions without reference to mixing- induced CP 
violation. A discrepancy between the constraints on 7 from charmless hadronic B decays 
and those from B and K mixing might arise from new physics contributions to either B and 
K mixing or the b —>■ s or b —>■ d penguins. 

Global analyses of charmless hadronic B decays in the framework of Q CD-improved 
factorization 0, |^ find a value for 7 of approximately 80°. However, these fits predict 
smaller branching fractions for B — > K*^tt^ and K *^7r^ than are observed experimentally, 
and removing these modes from the above analyses improves the fit quality. It has been 
suggested that the B — > K*7i modes may receive dynamical enhancements not accounted 
for in Refs. 0] and Our analysis focuses on B ^ K*^'k'^ with input from a modest 
number of other B — »• PV decaySjthus providing a complement to the global fits. 

Following the notation in Ref. [2] for SU(3) invariant amplitudes, we denote color- allowed 
tree amplitudes by t and gluonic penguins by p. Amplitudes for \^S\ = 1 transitions are 
primed, while those for AS* = transitions are unprimed. A subscript P or V indicates 
whether the spectator quark hadronizes into the pseudoscalar or vector meson, respectively. 
Since B K*^K^ is dominated by penguin annihilation and ly-exchange contributions 
or rescattering effects, and these decays have not been observed experimentally, we neglect 
such amplitudes in our analysis. The transition amplitude for B — > K*^7i^ is A{K*^7i'^) = 
— {p'p + t'p). The amplitudes t' and p' carry the CKM matrix elements V*fyus and V^^yts 
with weak phases 7 and tt, respectively. The amplitudes for the two charge states are given 
by 

A{K*+Tx^) = \p'p\-\t'p\e'^e'^ (1) 
A{K'"Tx+) = \p'p\-\t'p\e-'^e'\ (2) 

and we can express the CP- averaged amplitude as 

= \p'p\^ + \t'p\^ — 2 \pp\ \t'p\ cos7cos(5. (3) 

We identify squared amplitudes, \A\^ = A* A, with branching fractions, B, and we absorb 
all numerical factors, like Gp, mp, phase space integrals, decay constants, form factors, and 



A{K*+7r- 



A{K*-TT- 
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CKM matrix elements, into the definitions of tlie amplitudes. Most of the B branching 
fraction measurements in the literature are calculated assuming equal production of charged 
and neutral mesons. We correct these branching fractions by the ratio of B^B~ to B^B^ 
production rates, /+_//oo, as well as by the ratio of charged to neutral hfetimes, Tj^/tq. 
Because the constraints on 7 are constructed from ratios of branching fractions, we scale 
only the neutral B branching fractions by the product JF = ■ which is measured 
directly in Refs. and jij. Thus, cos 7 cos 5 can be expressed in terms of the CP-averaged 
branching fraction B{K* 



COS 7 cos 5 



\p'p 
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2 \p'p\ \t'p 
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For a given value of 5, 7 is determined to a twofold ambiguity. The rate difference between 
B^ K*~TT~^ and B^ K*~^7i~, which is proportional to sin 7 sin 5, provides an additional 
observable that allows us to disentangle 7 and 6: 
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These expressions for 7 and 5 are subject to a fourfold ambiguity: {7, 5} {5, 7}, {—7, —5}, 
or {-5, -7}. 

The charge-separated branching fractions B{K*^t[~) and B{K*^ -k^) appearing in Eqs. 
IHl can be determined directly from B{K*^tt^) and the CP asymmetry Acp{K*'^'k^). In 



addition, SU(3) symmetry relates the rate difference ^[K* tt"*") = B{K* 



to the corresponding AS' 



quantity, A(p-7r+) = B{B° 



p 7r^ 
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-77+) - B{K*+t:-) 
P+vr^) 0,113: 

(9) 



To attain greater precision on B{K*^it~) and B{K*^n^), we combine information on 
A(p"7r+) with the measurements of B{K*^7i'^) and AcpiK*"^?!-^) hsted in TableHl These in- 
puts are given relative weights that minimize the uncertainties on B{K*^7t~) and B{K*~tt~^), 
and we account for the correlation between CLEO's Acp{K*^tt^) and B{K*^tt^) measure- 
ments, which are made with the same dataset and technique. 

The BaBar analysis oi B ^ tt+tt^tt^ which determines the CP asymmetry and 



dilution parameters A, C, and AC defined in Ref. allows us to evaluate A(p vr 
-{A + C + AAC)-Bip 



We propagate the uncertainties on these parameters with their 



correlations jl2|J to obtain A(p 
in Eq.Elto be unity, we find B{K* 



n 



-(2.9±4.6)-10-6 
-) = (14.411^) -10- 



Thus, taking the form factor ratio 
' and B{K*-n+) = (18.7^) ■ IQ-^ 
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with a correlation coefficient of 0.53. The correlation coefficients between A(p~7r"'") and 

these two branching fractions are 0.45 for B{K*~^tt^) and —0.50 for B{K*~tt^). 

A second method of estimating 7 uses B — > K*^tt^ and B (pK^. The possibility 

of constraining 7 from these decays was ffist noticed by Gronau and Rosner 113^ , an d the 

concrete formulation of this method was subsequently put forth by Gronau |l5| . The 

SU(3) decomposition of the B (pK^ amplitude is A{(j)K^) = p'p — |-Pew where 
p/p 

EW 



denotes the electroweak penguin contribution. The weak phase of -P^w same as 



that of p'p, and its strong phase is expected to be the same as in t'p because of the similarity 
of their flavor topologies 3 E|- Thus, the ratio of the CP-averaged branching fractions 
for B K*^TT^ and B (j)K^ provides a measure of 7 up to a twofold ambiguity: 
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where r = \t'p/p'p\, and 



R 



\A{K*+'n-)\^ + \A{K*-'K+f 
\A{cl>K+)\' + \A{ct>K-)\' 



Both B <pK^ and B (f)^K^° receive the same SU(3) amplitude contributions 
can improve the statistical precision of Eq. ^2 by combining both channels: 
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so we 



(12) 



where 0"+ and ctq refer to the uncertainties on B{(t)K^) and B{(j)^K'^)T, respectively. To 
determine 7 with this method, the size of 5 must be known. It is believed, based on 
perturbative and statistical fl7^ calculations, that 0° < ML< 90°. In the simulation, we 
fix I Pew I to be | |pp|, as given by factorization calculations [3, IE 3, and we evaluate the 
dependence of our results on \P^/p'p\ and 5. 

In both of the above methods of constraining 7 (involving Eqs.EHHland Eg . ITIH) . numerical 
values of \t'p\ and \p'p\ are given by other B PV branching fractions [20]. The penguin 
amplitude is simply 

\p'p\ = ^JbCK^*^^). (13) 

The tree amplitude is taken from the AS = transition B p^vr^ and related to the 
I AS* I = 1 amplitude through SU(3)-breaking factors: 
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The experimentally measured Blp'^n^) represents a sum over B^ ^ p 
decays: 
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^13.9 ±2.7) • 10-*' and a 



Based on the experimental inputs in Table HI we find 
correlation coefficient between B{p^7r^)p and A(p~7r"'") of 0.05. 

Extracting \tp\ from B{p'^n^)p requires estimates of the magnitude and phase of pp. Its 
magnitude is obtained from the analogous \AS\ = 1 amplitude: 
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In the SU(3) limit, pp and tp have the same relative strong phase as that between p'p and 
t'p. Their relative weak phase, however, is 7 + /3, where 7 is unknown, a priori. Therefore, 
we must solve for cos 7 and \t'p\ simultaneously. 

Using CKM unitarity, the parameters |V^d/VJs| and j3 can be eliminated in favor of 
and 7 via the relations 
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By making these substitutions, we remove our dependence on sin 2(3 measurements involving 
b ^ c transitions and B^-B^ mixing, and we remain sensitive to new physics which may 
affect these processes and charmless b ^ u,d,s transitions differently. 

From the above unitarity relations and the CP-averaged branching fraction 



\ppf + \tpf + 2 \pp\ 



\tp\ cos(7 + (3) cos 6, 



(22) 



we find the following expression for \t'p\: 



where 
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Using Eq. Enito calculate \t'p\ from B K*^ti^, K*^'k^, and p^T[^ depends on a choice 
of 5 as well as knowledge of 7, and an iterative solution is required. The fixed strong phase 
appearing in Eq. |^ is distinct from the strong phase in the simulated quantities cos 7 cos 5 
(Eq. HJ and cos(7 ± 5) (Eqs. ElandlHl). To distinguish these two strong phases, we denote 
the one entering Eq. |^by 5^^. Below, we verify that the simulated values of cos 7 cos 5 and 
cos(7 ± 5) are insensitive to the choice of dt'^. In the second method of constraining cos 7, 
we simulate Eq. [TUlwith = 5. 

Experimental measurements of the following quantities are given as input to the simu- 
lation: fx*/ fp, \Vus\, \Vudl iKfel, \Vcbl -4cp(-ft'*^vr=F), the B p^ix^ parameters A, C, 



— \ Vus \ COS 7^ cos (5. 



(23) 
(24) 



and AC, and the CP- aver aged branching fractions for B 
and 



''K^'^. These parameters are simulated with Gaussian or bifurcated Gaussian (different 
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Parameter 




References 




Value 


T 




[7,8] 




1.11 ±0.07 


Ik* I fp 




[20] 




1.04 ± 0.02 
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[2JJ 




0.2196 ± 0.0026 
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[2JJ 




0.9734 ± 0.0008 
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(3.66 ±0.51) • 10-3 
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(4.07 ±0.10) • 10^2 






[22] 




n oc+0.33+0.10 
"^•^0-0.34-0.08 


A 




lUJ 




-0.18 ±0.08 ±0.03 


C 








0.36 ±0.18 ±0.04 


AC 




lUJ 




0.28l|j:J^ ± 0.04 


PAC 




[12] 




-0.080 


PAAC 




[12] 




-0.059 


PCAC 




fl2] 




0.176 




CLEO 


BaBar 


Belle 




V> { + =1= \ 

D(/3^7r+) 


27.6+^:^ ± 4.2 [23] 22.6 ± 1.8 ± 2.2 [11] 


20.8+^.^^3 5 [24] 


22.8 ± 2.5 


15\K vr j 


7.6l|g ± 1.6 [23] 


15.5 ± 1.81JJ [2^ 


1 n 0+4. 2+4.1 ro^;! 
19.3l3,9l7,i [2^ 


12.4 ± 2.5 




|l6l^ ±2 [22 
5.5l^:^ ±0.6 [30] 


lO.Ol^l ± 0.5 [31] 


20.3t^| ±4.4 [2^ 

iq n+3-9+2.0+6.9 rr,Qi 
J^'J-U_3. 6-1. 8-6.1 1^ 

10.7 ± 1.0+?-^ [32] 


}i6.4;tg 

9.6 ±0.8 


B{<i)''K^^) 


5.4l|^ ±0.7 [30] 


7.61};^ ± 0.5 [31] 


lo.oi}:?;?! [32] 


8.1 ± 1.1 



'8(p^7rT)p [11, 12, 23, 24] 13.9 ± 2.7 

A(p-7r+) [ll>J^>^,Ji] -2.9 ±4.6 

B{K*+Ti-) ^,XL2^^2^^M^2i^ZL2^^M i4.4tt^ 

B{K*-TT+) [11, VI , 2^, 22 , 2^, 24 , 2^, 2^, 2^ 18.71^6 



TABLE I: Input parameters used to constrain 7. Branching fractions and partial rate differences 
are given in units of 10"^. Except for the last two entries, branching fractions are averaged over 
charge conjugate states. 



widths above and below the peak) distributions, and their values are summarized in Ta- 
ble H] The input that contributes the largest uncertainty to 7 is the B JC*^7r^ branching 
fraction. 

For the five branching fractions, we combine all publicly presented measurements, with 
statistical and systematic errors added in quadrature. Where possible, the contribution 
from /_|_-//oo to the systematic error has been removed, since it is included coherently in 
the simulation. We neglect all other correlations among the systematic errors. 

Figure HI shows the simulated distribution of cos 7 cos 5 from Eq. HI with = 0°. Fitting 
this distribution to a bifurcated Gaussian yields the measurement cos 7 cos 5 = — O.GS+olg, 
which suggests constructive interference between t'p and p'p. The 90%, 95%, and 99% C.L. 
upper limits on cos 7 cos 5\5 , =0° are 0.16, 0.42, and 0.94, respectively. Based on the smallness 

of direct CP asymmetries in i? ^ Ktt, one can infer a strong phase between tree and penguin 
amplitudes in these decays of (8 ib 10)° js^J- If the strong phases in i? ^ PV decays are as 
small as in two-pseudoscalar {PP) final states, then our analysis weakly favors cos 7 < 0. 
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FIG. 1: Simulated distributions of cos7cos5 from Eq.3][(a), solid squares], cos(7 + 5) from Eq.|Sl 
[(b), solid squares], and cos(7 — S) from Eq. El [(b), open circles] using B K*^-7t^, as well as 
COS7 from Eq. El[(a), open circles] using B K*^tt^ and B (j^K^ with \P^^/p'p\ = 0.5, all 
with S^i^ = 0°. Overlaid on the histograms are the fits to bifurcated Gaussians. The dashed lines 
demarcate the physical region. 



The variation of cos 7 cos 5 with cos5f' is roughly linear, with a slope of = 0.11. 

Also shown in Figure Q is the distribution of cos 7 from Eq. with 5 = 5t' = 0° and 
\Pl[^/p'p\ = 0.5. Here, we obtain cos 71^=0° = -O.bOtoM and 90%, 95%, and 99% C.L. 
upper limits on cos 7(5=0° of 0.23, 0.44, and 0.89. 

From Figure m we also find cos(7 + 6) and cos(7 — 6) from Eqs. ElandlHlto be — 0.39lo;g3 
and -0.99tlit, respectively, with a correlation coefficient of 0.61. Considering only the 47% 
of trials where both quantities acquire physical values, the distributions of the weak and 
strong phases imply 7 = (118^30)° and 5 = (—13 ± 17)°, with a correlation coefficient of 
7 ■ 10~^. Because of the fourfold ambiguity of the 7/5 system, we fix 6t'^ to 0° rather than 
equating it to the simulated value of 6. The variations of cos(7 + 6) and cos(7 — 6) with 
cos 5+' are given by = q qq and = 0.13. The values of 7 and 5 both change 

'•p ° dcosOj/ a cos 0^1 ' ° 

by less than 2° between 6t> = 0° and 6t> = 80°. 

Figure El shows the dependence of C0S7 from Eq. (TUlon 6 = 6t' , with \Pew/p'p\ — 
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FIG. 2: Peak values (a) of and upper limits (b) on cos 7 from Eq.^Jas a function of 5 = (5j^, using 

B — > K*^iT^ and B — > (p^'K^ with \P-ew/p'p\ ~ ^-^^ '^^^ asymmetric errors on the peak values 
give the bifurcated Gaussian widths of the simulated distributions. The dashed lines demarcate 
the physical region. 

The peak values are plotted with asymmetric error bars representing the widths of the 
bifurcated Gaussian distributions. By demanding that cos 7 peak in the physical region, 
one can infer that \6\ < 41°. The variation of C0S7 with |-Pew/PpI Hnear, with a slope 
I '^^p^/' , I = 0.28 — 1.51 cos 6t' . Incorporating the B — > (j)^K^ decays in the measurement of 7 

results in greater precision than using B K*^7t^ alone, but the theoretical uncertainties 
incurred are also larger. 

Using the simulation of Eq. 0J we also determine the ratio r = 0.301q;o5 at 6t'^ = 0° with a 
St'^ dependence given by r = 0.25 + 0.09 cos 61'^ — 0.04 cos^ 6t'^ . The inverse ratio for AS* = 
decays, \pp/tp\ = ^ \Vtd/Vts\, is found to be 0.43-0.49 cos 6t' +0.26 cos 6t' , which 

takes the value 0.20lo:o2 at St'^ = 0°. 

The widths of the generated distributions presented above are dominated by experimental 
uncertainties on the input branching fractions, Acp{K*'^7i'^), A, C, and AC. We study the 
improvement in the resolutions of cos7cos5, cos(7±(5), and cos7|5=o°) collectively denoted 
by o"cos7, as these measurement uncertainties are reduced while maintaining the central values 



8 



at their current positions, with St'^ = 0° and |-Pew/PpI — It is found that (3"cos7 scales 
with the size of the experimental uncertainties until the latter reach 10% of their current 
values, where the resolution of 7 is 0(10°). At this point, (Tcos7 begins to be dominated by 
the uncertainty on JF, and only by lowering ajr can (3"cos7 be reduced any further. 

We have modeled SU(3) symmetry breaking effects in ratios of AS" = to \^S\ = 1 
amplitudes with the purely real ratio of decay constants fx*/ fp- Repeating the simulation 
without SU(3) breaking (z.e., with fx*/ fp = 1) results in changes to cos 7 cos 5, cos(7 ± 6), 
and cos7|5=o° of 0.05 or smaller. Recent studies based on QCD-improved factorization [lO|, 
S^l have suggested that SU(3) breaking could be as large as 30% and that the amplitude 



ratios may possess a small complex phase. To probe the impact of such effects, we reinterpret 
fx*/ fp as a phenomenological parameter and scale it by ±30% of the value given in Table HI 
neglecting any possible complex phases. We find shifts of cos 7 cos 5, tois in cos(7 + 

5), -o.il cos(7 — 6), and taso in cos7|5=o°- Thus, in this conservative estimate, SU(3) 
breaking effects are roughly 15%-70% of the current experimental uncertainties. To obtain 
meaningful constraints on 7, future experimental advances must be accompanied by an 
improved understanding of SU(3) breaking. 

In conclusion, we have formed constraints on 7 as a function of 6 and \Pew/Pp\ using 
branching fractions of and CP asymmetries in B ^ PV decays. At present, experimental 
uncertainties overwhelm the theoretical uncertainties arising from the model dependence 
of \Vub\ and \Vcb\, but they are the same order of magnitude as the uncertainties in SU(3) 
symmetry breaking. For strong phases of (9(10°) or smaller, our analysis favors cos 7 < 0, 
which agrees with indications from B — * PP decays [s^, 11^ 3^. However, the current 
experimental precision does not yet permit a stringent comparision with fits reliant upon B 
and K mixing. 
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